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By Richard S. Cesaro, Robert J. Koenig 
and George J. Pack 


SUMMARY 

Various theoretical and practical aspects of a pressure - 
sensitive system for the measurement of gas temperatures in gas- 
turbine -engine combustion chambers are analyzed. An experimental 
setup was operated under controlled conditions of temperature , 
weight flow, and wall temperature. The gas temperature after com- 
bustion was obtained by application of a relation equating thermo- 
dynamic conditions before and after combustion. This method of gas- 
temperature sensing appears practical for high -temperature appli- 
cations. The accuracy with which gas temperature may be determined 
by this method is within ±2 percent. Measurements made by the 
thermodynamic method were used as a temperature standard for com- 
parison of temperature data obtained from conventional thermocouple 
probes . 


INTRODUCTION 

With the advent of gas-turbine engines, the controls designer 
has been confronted with the serious problem of operating engines 
not only at high temperatures for maximum efficiency, but also at 
temperatures low enough to be within the safe operating range for 
the materials used. 

It therefore becomes imperative to incorporated temperature - 
limiting controls to prevent engine failure by overheating. 
Temperature can be used both as a primary control parameter and as a 
temperature -sensing parameter. Any means of sensing gas temperature 
for application to controls must meet the following requirements : 

(1) have an output that is a simple function of true gas temperature 
under various conditions of weight flow, temperature of surrounding 
materials, density, composition of gas, and temperature level; 

(2) provide an output that is easily incorporated in a control; 

(3) have a rapid response to transient conditions; (4) have an 
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extended temperature range; (5) have a long life; and (6) be 
independent of material characteristics subject to change with 
time and uae. 

The most commonly used methods for measuring temperature 
employ thermocouples or resistance -wire thermometers, both of which 
have the advantages of simplicity and durability. These devices, 
however, do not provide sufficiently ideal outputs when consider- 
ation is given to the inherent errors resulting from radiation, 
heat conduction, and velocity losses that occur at elevated gas 
temperatures . In addition, chemical action on instrument materials 
can result in a change of calibration as well as in a reduction in 
life, which is a particularly difficult factor to contend with 
because the mass of the unit must be large in order to lengthen life; 
whereas the requirement of rapid response to a transient demands that 
the mass be small . At temperatures below 1000° F these errors can 
be reasonably neglected in control work;. At elevated temperatures, 
however, such errors cannot be neglected, particularly radiation 
errors, because these errors appear as the difference of the fourth- 
power functions of the two temperatures involved. 

As a possible means of circumventing present difficulties, a 
theoretical analysis was conducted at the NACA Lewis laboratory 
using the thermodynamic properties of gases as basic parameters. An 
equation has been developed based on these properties (reference 1). 
This equation, based upon the expression for weight flow, correlates 
gas temperatures with measured pressures before and after combus- 
tion. Such pressure measurements are independent of radiation or 
heat -conduction effects, and therefore permit an evaluation of the 
true gas temperature (within the accuracy of pressure -measuring 
instrumentation) as long as the ideal gas laws apply. 

Various theoretical and practical aspects of the pressure - 
sensitive system for gas -temperature measurement are discussed herein. 
The system was applied to an experimental burner setup operating 
under controlled conditions of temperature level, weight flow, and 
test -section -wall temperature. Several methods of measuring weight 
flow were analyzed to determine their applicability to the system. 
Particular emphasis is given to a discussion of the pitot -static 
method of determining weight flow by analyzing the velocity profile 
across a test section at elevated temperatures . 

For comparison purposes, the pressure -sensitive method (refer- 
ence 1) is considered as a reference system for sensing temperature. 
Data obtained with several conventional thermocouple probes were 
compared with data from this reference system over a range of 
temperature levels, weight flows, and duct-wall temperatures. 
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ANALYSIS 

Various types of fluid meter, including orifices, nozzles, 
venturis, and pitot tubes, or combinations of such units, may be 
employed to evaluate thermodynamic relations of reference 1 in 
terms of gas temperature. Experimental data reported herein were 
obtained using an orifice before combustion (station l) and a 
pitot-static tube after combustion (station 2). The general 
weight-flow equation, which applies for all fluid meters, is 


W = pAV 


( 1 ) 


(Definitions of symbols used in this section are given in 
appendix A. ) 


Orifice Equation 

The specific equation for flow through an orifice, developed 
from equation (l), is 


W x = A 1 K L E 1 Y 1 '\/p 1 AP 1 2g (2) 

The measurements required by equation (2) for the determination of 
weight flow at station 1 can be readily and accurately obtained 
using an orifice installation, as recommended in reference 2. 


Pitot -Tube Equation 

Discussion of pitot -tube method . - A pitot -static tube was 
used at station 2 because of the convenience of this method when 
dealing with high -temperature gases, inasmuch as it offers negligible 
resistance to gas flow, has a rapid response (reference 3) , and is 
simple in construction. The specific equation for weight flow when 
applying the pitot -static tube for flow measurement at station 2, as 
developed from equation (l), is 

W 2 = ^ 2 ^ 2 ® 2 ^ 2g 2 P 2 AP 2 ^ 

When using this method it is important to realize that the 
measured value of AP (the total pressure minus the static pressure) 
is that of the specific point at which the pressure -sensing probe is 
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located. In general, the flow or velocity profile across a section 
is not flat, therefore AP is not constant across that section. In 
order to make the system applicable to the temperature equation, it 
is necessary that a single -point indication he proportional to the 
total weight flow through the section. It therefore becomes essential 
to determine the relation that exists between a measured value of 
AP, which is quite readily obtained, and the effective value of AP, 
which is a function of the average velocity of the gas . If it is 
assumed that the effects of changes in boundary layer on weight flow 
over the range of temperature operation considered remain essentially 
constant, the ratio of the effective AP to a measured AP can be 
evaluated. 

In order to successfully apply equation (3) for a given 

apparatus, whether it be an aircraft power plant or a simple straight- 

through pipe, the relation between a measured AP and the effective 

AP must be consistent (AP „„/AP m = C, a constant). If this 

six ms a. s 

consistency does not exist, the method presented herein for gas- 
temperature evaluation cannot be applied. It should be noted, however, 
that a consistent ratio of ^P e ff/^meas *> e considered as such 

only in regard to the individual application and the accuracy of gas- 
temperature evaluation desired. For example, in one installation 
an evaluation of gas temperature to within ±10 percent may be suf- 
fient; the ratio AP e ff/AP meaa therefore may vary ±10 percent and 

be considered consistent. In another installation, where a gas- 

temperature -evaluation accuracy of ±2 percent is desired, a ±3 percent 

variation in AP „,,/AP could not be considered consistent. For 

sir ms S.S 

ths data pressntsd, ths svaluation of A^eff/^meas was mads ovsr 

an extended temperature range for fully developed turbulent flow in 
a pipe; furthermore, AP was measured at the center of the gas stream 

SO that hPjjjgg^g = hPjjjgjj-. 

Evaluation of AP pff /AP^^ . - If a fully developed laminar- 

flow velocity profile is realized in the section, it is possible to 
determine the ratio AP p ,- f .. P /AP msT from theoretical considerations. 

Because experimental data were obtained for turbulent -flow conditions 
only, experimental methods had to be employed in order to evaluate 
the ratio. 

The velocity profile in a duct cannot be defined by a simple 
mathematical relation when turbulent -flow conditions are obtained. 

The profile is more nearly uniform under these conditions than for 
laminar flow. Figure 1 shows a typical curve of the velocity ratio 
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v/v ’ as a function of Eeynolds number. These data were recorded 

mai 


by Stanton and Pannell in 1914 and were obtained from reference 4. 
The curve shows that in the laminar region A the velocity ratio 
is of the order 0.5. A very rapid change from 0.5 to 0.72 is noted 
in the transition region B. In the fully developed turbulent 
region C , the velocity ratio is nearly constant at a value of 
approximately 0.78. It should be noted, however, that these 
numerical values probably would not apply to flow conditions that 
exist in turbojet engines. The important point is that a consistant 
velocity ratio must exist before a single pitot tube can be used to 
determine AP ef . f for any application. 

A direct method of evaluating AP eff /AP max is to obtain 

actual AP profile data by making a traverse of the section using 
a pitot tube . Particular attention must be given to recording 
accurately data obtained near the pipe walls, because a steep 
gradient exists in that region. 

Inasmuch as 


where V, the average gas velocity, is a function of ^AP ef . f 

when p is a function of p/BT (appendix B), and A is a function 
of radius squared, a value of AP^^. can be obtained from consider- 
ations of AP (measured). In order to obtain this value, it is 
necessary to plot the square root of the measured pressures against 
the square of the radius . Figure 2 shows that /\/AP ^ is propor- 
tional to A', the area under the curve, and f\J A P mn x is propor- 
tional to area A 1 plus A", the hatched area. Therefore 


The ratio of AP g ff /AP max can be determined quite readily 
from an analysis of measured pressures. 


W = p AY 


( 1 ) 



or 
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The average AP, calculated from a traverse of the section 
with a pressure probe, does not have the same significance as the 
effective AP based on the considerations of an average velocity 
across that section. 

Another method of evaluating the ratio AP e ff/AP max can be 

used if the gas temperature after combustion (station 2) is known 
for at least part of the temperature range. For this case, the 
weight flow at stations 1 and 2 can be equated and solved for 
AP e pf; the ratio can then be found by using a measured value of 

AP ma;x . In order to equate the weight flow at stations 1 and 2, it 
is necessary to account for the addition of fuel: 


W 2 = W x + f = 



(5) 


The weight flow at station 2 (equation (3)) is equal to the 


weight flow at station 1 (equation (2)) times 



When a 


thermocouple is used to determine gas temperature at station 2, 
as it was for this experimental analysis, errors can be expected at 
temperatures exceeding 1000° F, at which temperature radiation errors 
of thermocouples become prominent. The weight-flow evaluation 
method was used as a check on the pressure -survey method only in the 
range in which errors of thermocouples could be neglected. A 
theoretical discussion of thermocouples and their inherent errors is 
given in appendix C . 


Temperature Equation 

General equation . - The orifice and pitot -static -tube equations 
for weight flow may be equated: 


E 2 A 2 C f l ' 2 / \l 2 8 p 2 AP 2 = A 1 K 1 E 1 Y 1 V p l AP l 2 3 [ 1 + ~''j 

When the density p is replaced by its equivalent JL (appendix B), 

KT 

equation (6) can be solved for Tg . 
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(7) 


Equation (7) represents the complete temperature equation with 
which experimental data discussed herein will he concerned. 

Simplified temperature equation . - In reference 1 it is shown 
that, based on analytical considerations, the term 


of equation (7) is very nearly constant. Considering the term as a 
constant, the temperature equation becomes 


For the results presented herein, each factor of the coefficient N 
has been taken into consideration and the analytical conclusions 
substantiated . 


A schematic representation of the apparatus used to obtain 
the data presented is given in figure 3. Air was drawn from the 
atmosphere and compressed to a pressure of approximately 3.5 pounds 
per square inch gage by two commerical blowers connected in series, 
then ducted through the combustion chambers and test sections. The 
maximum air flow available was approximately 6.5 pounds per second 
at a pressure of 2.5 pounds per square inch gage in the test 
sections. An orifice designed and installed according to A.S.M.E. 
specifications (reference 2) was located at station 1, downstream 
of the blowers, for the determination of air flows. Located down- 
stream of station 1 were four combustion chambers capable of heat- 
ing the maximum air flow to a temperature of 2500° F. A flame 
shield and gas mixer, installed to shield the thermocouples in 
test sections from radiation from the flame front and also to 
effect uniform gas mixing, was located immediately downstream of 
the combustion chambers. The gas leaving the flame shield and gas 
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mixer was ducted through a 10-inch-diameter Inconel pipe for 
approximately 15 feet, whereupon it entered test sections 1, 2, 
and 3, also 10 inches in diameter and arranged in series, and 
was finally discharged to the atmosphere . 

The instrumentation at test section 1 (fig. 4), designated 
station 2, consisted of: a six -tube static -pressure survey rake; a 
six-thermocouple survey rake incorporating single -shielded thermo- 
couples; three thermocouples embedded in the wall of the test 
section; and a total -pres sure rake. A total -head tube of 0.030- 
inch outside diameter, which could be positioned along the pipe 
diameter during a run for the purpose of obtaining a complete 
total -pressure -profile survey was also used. 

Test section 2, illustrated in figure 5, contained eight 
thermocouples having bead diameters of 3/8, l/4, 3/l6, 5/32, l/8, 
3/32, l/l6, and l/32 inch, which were arranged in a ring concentric 
with the pipe. Three additional thermocouples were embedded in the 
wall of the test section. 

Test section 3 (fig. 6) contained the following instrumentation: 
a thermocouple having an electrically heated shield; a commerical- 
type multiple -shield thermocouple; a gold -encased single -shield 
Bureau of Standards thermocouple (reference 5); three thermocouples 
embedded in the wall of the test section; and a pitot-static pressure 
tube in the gas stream. A cutaway drawing of the heated -shield 
thermocouple is presented in figure 7. During steady-state operation 
the electric -heater coil supplied heat to the shield surrounding the 
thermocouple located in the gas stream until the indicated wall 
temperature of the shield equalled indicated gas temperature . When 
this equality existed, the gas temperature was recorded under the 
assumption that conduction and radiation losses were then minimized. 

Operating conditions were established by maintaining the weight 
flow at a constant level and varying the fuel flow to the combustion 
chambers, whereby the time required for changes from one steady- 
state temperature level to another was held to a minimum. In order 
to present the effects of weight-flow variations on the operation of 
the various temperature -sensing probes, three weight flows were inves 
tigated over a temperature range from approximately 900° to 2200° F. 
The flows selected for these tests were approximately 4.5, 5, and 
6.5 pounds per second. A plot of weight-flow variations at station 2 
is presented in figure 8. For the temperature range considered, the 
gas velocities in the test sections varied between 300 and 800 feet 
per second with the various weight flows. The Mach number, however, 
varied only slightly for each setting, as can be seen from figure 9. 
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The effects of conduction and radiation losses on the thermo- 
couple indications of gas temperature at constant ref erence -gas - 
temperature settings were investigated by controlling the test- 
section-wall temperature with cooling-water sprays acting on the 
outside walls of the test sections. For runs in which water cooling 
was not used, there was no variation in indicated wall temperature 
at any given reference gas temperature during variations of weight 
flow. This effect is illustrated in figure 10. 


RESULTS AND DISCUSSION 
Evaluation of Constant C 

Profile survey . - The results of pressure -profile surveys 
taken at station 2 for reference gas temperatures of 550°, 1500°, 
and 2300° R are presented in figure 11. A plot of C for all 
operating gas temperatures is given in figure 12. Any differences 
in boundary layer caused hy differences in gas temperature (ref- 
erence 6) are small enough to cause a variation in C of less than 
0.5 percent. On the basis of these results, a value for C of 0.87 
was selected and used in the evaluation of the temperature equation 
(equation (7)) as applied to data presented herein. 

Weight -flow equations . - The evaluation of C obtained by 
equating the weight flow at stations 1 and 2, using thermocouple 
data taken at these stations, is presented in figure 13, As 
mentioned in the section ANALYSIS and discussed in appendix C, inherent 
errors of the thermocouple should be considered at temperatures of 
approximately 1000° F (1460° R) and above, which is shown by the data 
presented in figure 13. A comparison between the weight-flow and the 
pressure -survey evaluations of C is shown in figure 14. 

In view of the results obtained by the two methods of evalua- 
tion of C, and the simplicity of weight-flow method over the pro- 
file survey, the weight-flow method should be satisfactory for use 
in determining C in practical application of pressure -sensitive 
temperature -sensing systems, providing this determination is made 
at temperatures below 1000° F. 


Pressure -Sensitive System as Temperature Standard 

On the basis of the fundamental relations used in deriving the 
temperature equation (equation (7)) and the experimental evaluation 
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of the constant C, the pressure -sensitive system is hereinafter 
considered a temperature standard. This standard, "based upon the 
application of equation (7), introduces only errors of instrument 
reading; the temperature -sensing accuracy is therefore considered 
to "be within ± 2 percent . 

Weighing is the fundamental method of fluid measurement and is 
the technique used hy the A.S.M.E. to establish a flow-measurement 
standard for fluid meters, such as an orifice plate, pitot -static 
tube, nozzle, venturi, and so forth (reference 2). The A.S.M.E. 
has established flow -measurement accuracy of approximately ±1 per- 
cent for these various flow-measurement units. Installation and 
instrumentation requirements, as established by the A.S.M.E. in 
order to obtain this accuracy, were used in obtaining the data for 
the evaluation of the flow equations presented. Inasmuch as the 
maximum error in weight-flow measurement is approximately ±1 percent, 
in effect, the maximum deviation of any one parameter in equation (3) 
is limited to ±1 percent, and as weight -flow is evaluated at two 
separate stations, the maximum error deviation between the two 
stations is ±2 percent. Because density, and therefore temperature, 
is one of the parameters in equation (3), the greatest deviation of 
temperature is ±2 percent from a true absolute level. The state of 
development of pressure -measuring techniques is such as to insure a 
rapid response rate for the pressure -sensitive system (reference 3) . 
The error involved in the use of the simplified form of the tempera- 
ture equation (equation (8)) is shown to be quite small in figure 15, 
where values of the dimensionless coefficient N, plotted against gas 
temperature for all operating conditions, fall within a band of ±1.5 
percent. Equation (8) was not used in connection with the temperature 
standard . 


COMPARISON OF THERMOCOUPLES WITH 
TEMPERATURE STANDARD 

No attempt is made to evaluate the various thermocouples in 
this discussion. A general explanation of the performance of the 
various thermocouples is given in the theoretical, discussion of 
thermocouple operation in appendix C . 


Radial Temperature Profile 

A typical plot of gas -temperature profile, as measured by the 
single -shielded -thermocouple survey rake in test section 1, is com- 
pared to the temperature standard in figure 16 for reference gas 
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temperatures of approximately 1450° and 2000° R. Because the gas 
temperature is uniform across the section, a comparison can he made 
between any thermocouple -temperature indication and the reference 
gas temperature given by the temperature standard or reference 
system. 


Performance of Thermocouples 

A performance plot is presented in figure 17, showing tempera- 
ture deviations from the reference temperature of (l) a commercial 
multiple shield, (2) a single-shield survey rake across the pipe 
diameter (average of all temperatures), (3) a Bureau of Standards 
gold shield, (4) a KACA heated shield, and (5) a 3/8-inch-diameter 
junction thermocouple for weight flows of 4.6, 5.3, and 6.6 pounds 
per second. Bata for the Bureau of Standards thermocouple were not 
obtainable for the mass flow of 6.6 pounds per second because the 
thermocouple burned out during the performance run. 

A plot showing the percentage deviation from the calculated 
reference temperature for four of the five types of thermocouple is 
given in figure 18 for weight flows of 4.6, and 6.6 pounds per 
second. Prom the slope of the curves, it is apparent that increas- 
ing reference gas temperature considerably increases the percentage 
of error. 


Influence of Thermocouple -Bead Diameter 

Indications of gas temperature obtained from thermocouples of 
several bead diameters are compared with the reference temperature at 
2525° R in figure 19. The data were extrapolated in order to obtain 
an intersection with the reference temperature at a thermocouple - 
bead diameter of zero. Without a reference temperature to indicate 
the true end condition of gas temperature at a thermocouple -bead 
diameter of zero, extrapolation of the data curves could be greatly 
in error. Also, the differences in temperature indication of the 
same thermocouple operating at a constant reference gas temperature 
but at various weight flows are quite apparent from the data. An 
equation representing the data curve would be impractical when added 
consideration is given to such factors as varying wall temperatures, 
Mach number, and reference levels . 


Effect of Wall Temperature on Thermocouple Indications 

The effect on thermocouple indication of changes in conduction 
and radiation losses at various reference gas temperatures were 
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investigated by controlling test-section -wall temperatures with 
cooling -water sprays acting on the outside of the walls . A plot 
showing the deviation of several thermocouples from the tempera- 
ture reference for a wall temperature of 710° E is given in 
figure 20. (The effects on thermocouple indications of changes 
in conduction and radiation losses may "be seen when comparing 
figs. 20 and 17(c).) The percentage deviation is given in fig- 
ures 21(a) and 21(b) for weight flows of 4.6 and 6.6 pounds per 
second, respectively. The slope of the curves of figure 19 shows 
that large errors can be encountered at the higher reference gas 
temperatures . A tabulation of results giving the performance of 
the temperature -sensing probes investigated at reference gas 
temperatures of 2000° and 2400° E and various weight-flow con- 
ditions, is given in table I. 


CONCLUSIONS 

Various theoretical and practical aspects of the pressure - 
sensitive systems for the measurement of gas temperatures have been 
discussed. The gas temperature after combustion was obtained by 
application of a relation equating thermodynamic gas conditions 
before and after combustion. The following conclusions may be 
drawn from this application: 

1. Determination of gas temperatures after combustion from 
measurements of gas temperatures before combustion and gas pressures 
before and after combustion appears practical for high-temperature 
applications. The temperature -sensing accuracy of this method is 
within ±2 percent . 

2 . Temperature changes of the combustion gases are accompanied 
by pressure changes that are in effect instantaneous, are unaffected 
by errors associated with conventional thermocouples, and can be 
utilized in a control system. 

3. By using the thermodynamic method as a temperature standard 
for comparison of thermocouples, it is determined that conventional 
thermocouples can be appreciably in error in their temperature indi- 
cations at gas temperatures above approximately 1500° E. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, July 29, 1949. 
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APPENDIX A 
SYMBOLS 


The following symbols are used in the section ANALYSIS and 
in appendix B: 

A area, sq ft 

C c ons t ant , AP eff / AP mnv 


E 

f 

g 

K 

N 


area multiplier for thermal expansion 
fuel flow, lb/sec 

acceleration due to gravity, ft/sec^ 
flow coefficient of orifice plate 


constant. 



P total pressure, lb/sq ft absolute 

p static pressure, lb/sq ft absolute 

E gas constant 

r ratio of static to total pressure, p/p 

T total temperature, °E 

t static temperature, °K 

V average velocity, ft/sec 

W weight flow, lb/sec 

Y ratio of flow coefficient of gas to that for liquid at same 

Eeynolds number (reference 2) 

y ratio of specific heats at constant pressure and constant 

volume 

AP total pressure minus static pressure, P - p, lb/sq ft 
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AP e f f pressure of average velocity, l"b/ eq ft 
frE’max pressure of maximum velocity, lb/sq ft 
AP measured value of total minus static pressure, lb/sq ft 
p density, p/PT, lb/cu ft 

cp conversion factor of hydraulic equation to compressible - 

flow equation 

cp' particular value of cp when p = p/ET (See equation B.) 

Subscripts : 
max maximum 

1 station 1 (before combustion) 

2 station 2 (after combustion) 


1222 


NACA TN 2043 


15 


APPENDIX B 


COMPKESSIBLE -PLOW EQUATIONS 


The hydraulic equation for incompressible flov may he 
multiplied by an appropriate conversion factor cp to obtain the 
exact equation for compressible flow. The expression for the 
conversion factor cp may be derived from the compressible -flow 
equation by factoring out the hydraulic equation so that the 
remaining factor is the expression for the conversion factor cp . 

Bernoulli's theorem for compressible flow may be written 
as 




The weight flow is 


W = Ap 0 V 


(Bl) 


(B2) 


Substitution of equation (Bl) in equation (B2) and replacement of 
stagnation density p t with the equivalent adiabatic relation 



The free -stream density p Q may be replaced by its equivalent 
and the equation simplified: 
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This expression for compressible flow may be written as 


W = Acp^2gp (p - p) 


(B4) 


where 



1 ( 7 \ ]£. 

(p - p) \y - 1/ Rt 




1 

2 


(B5) 


This expression for cp involves the density p, which also appears 
in equation (B4) and may he arbitrarily selected as a ratio involv- 
ing total pressure P or static pressure p divided by total 
temperature T or static temperature t . For the case in which 
the density p is selected as p/RT, the conversion factor cp is 
designated cp' and equation (B5) is simplified as follows: 
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The adiabatic relation of the temperatures is 


T 

t 



(B7) 


This relation may be substituted into equation (B6) to obtain 



If the pressure ratio p/P is set equal to r this equation 
becomes , \ 


(cp') 2 - 


1-7 
7 , 


1-7 

,T 7 -1, 


? - 1J (7-D 


or 


fo ') 2 



(1 - r) (r - 1) 


(B8) 


Evaluation of the conversion factor cp* for various pressure 
ratios r show the error that may be expected from neglecting 
cp in the hydraulic equation (B4) where the density p is p ,/PT . 
The greatest deviation of the conversion factor cp' from 1 occurs 
at the critical pressure ratio, at which cp' is approximately 
0.945 for the ratio of specific heats 7 equal to 1.3 for air at 
a temperature of 3000° R. 
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APPENDIX C 


CHARACTERISTICS OF THERMOCOUPLES 

A "brief summary of the significance of thermocouple indications 
is required in order to explain and to compare temperature data 
obtained from thermocouples with data obtained from the reference 
system . 

When dealing with thermocouples it must be noted that the 
indicated temperature is that of the junction itself and not nec- 
essarily that of the medium in which the junction is immersed. A 
thermocouple probe inserted in a hot gas stream indicates the 
temperature attained by the thermal element when a condition of 
thermal equilibrium with its environment is reached. The indicated 
temperature is therefore a function of several modes of heat 
transfer, including convection, conduction, and radiation. 

Velocity effect . - Actual total temperature is also a function 
of gas -stream velocity. This relation can be expressed as 


T = T 0 + T a 


where 


■o 


T„ 


free -stream temperature 


temperature due to velocity, 


V 2 

2gJC p 


J mechanical equivalent of heat 


specific heat at constant pressure 


Standard probes are designed to partly stagnate the gas stream 
around the junction and thereby recover a part of the temperature due 
to velocity. Because probes do not recover all this temperature, 
they must be calibrated for this effect so that the corrected 
expression is 



where r is the recovery factor. 
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Inasmuch as the probe temperature is a function of "heat in" 
minus heat out, these terms must be evaluated and analyzed to 
determine their effect on the final indicated temperature . 

Heat in is basically a function of heat transfer by convection 
from the gas to the solid material of the thermocouple, which can 
be expressed as 


q c = Sh(T - T p ) 

where 

q„ "heat in" 

S surface ar6a in contact with gas 

h heat -transfer coefficient 

T temperature indicated by probe 

XT 

The coefficient of heat transfer is an involved function of 
probe area and weight flow across the thermocouple as well as the 
viscosity and the composition of the gas . In order for the probe 
equilibrium temperature to approach total temperature, the difference 
(T - Tp) must be small. It therefore becomes evident that in order 

to insure sufficient heat transfer by convection, the weight flow 
over the junction must be high. 

Heat out (heat lost) is due to: (l) conduction of heat away 
from the junction through the leads of the thermocouple and, 

(2) radiation of heat to cooler surfaces. (It should be noted 
that if the walls of the duct containing the gas are at a higher 
temperature than the gas, then all heat flow is reversed.) 

Conduction errors in indications are given by 



where 

T^ temperature at distance L from junction 
L some distance from junction along wires 
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c circumference of wires 

k cross-sectional thermal conductivity of thermocouple wires 

A cross-sectional area of wires 

D diameter of wires 

From this relation it follows that in order to keep the 
conduction error small the l/d ratio and h must he large, 
whereas k must he small . 

Heat loss due to radiation may he one of the most serious 
errors encountered in using thermocouples for sensing elevated gas 
temperatures. This fact is readily appreciated w"hen an expression 
governing radiation error is analyzed. 

The heat loss due to radiation may he expressed as 
H r = e 8 S (T p 2 * 4 - T w 4 ) 

where 

e emissivity of radiating hody 

8 Stef an -Bolt zman constant 

T wall temperature of radiating hody 

w 

It should he noted that for a given installation of a thermo- 
couple prohe, the only variables concerned with the radiation error 
are the temperature of the junction and the temperature of the wall 
or radiating hody. Slight changes in temperature appear as the 
fourth power, making the difference and the resulting radiation 
error quite appreciable . 
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TABLE I 


PERFORMANCE OF TEMPERATURE - SENSING PROBES 



Weight flow 
(lb/sec) 

Mach 

Number 

Velocity 

(ft/sec) 

Calcu- 

lated 

average 

gas 

tempera- 

ture 

(°R) 

Deviation of thermocouple from calculated 

gas temperature, percent 

Nonshielded 

thermocouple 

Single -shielded 
thermocouple 

Gold -shielded 
thermocouple 

Multiple -shielded 
thermocouple 

Without 

cooling 

With vail 
cooling 

Without 

cooling 

With vail 
cooling 

Without 

cooling 

With vail 
cooling 

Without 

cooling 

With vail 
cooling 

4.60 

0.196 

425 

2000 

7.50 

7.70 

5.50 

6.10 

4.20 


3.40 

3.10 

6.55 

.280 

610 

2000 

6.65 

8.50 

4.60 

7.30 




3.20 

4.70 

4.50 

.214 

510 

2400 

10.70 

13.00 

8.20 

10.60 

6.30 



5.00 

6.60 

6.35 

.300 

715 

2400 

8.75 

13.00 

5.70 

11.00 

— 



3.50 

7.00 
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Velocity ratio, V/V, 


±.UL.Lu 


I 1.00 

.80 


.60 

.50 

.40 

.30 






A Laminar-flow range 
B Transition region 
C Turbulent -flow range 
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Reynolds number 


Figure 1. - Effect of Reynolds number on ratio of average 
velocity V/V ma for isothermal flow in smooth pipes, 
from reference 4.) 


to maximum 
(Data obtained 
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Total minus static pressure, yS"P 
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Figure 2, - Effect of Reynolds number on flow profile for turbulent flow 
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Test section 



Figure 3. - Schematic representation of research apparatus. 
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Figure 4. - Test section 1; pitot-static and thermocouple survey rake and instrumentation, 
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Figure 5. - Test section 2; varying thermocouple-bead diameter 
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Figure 6. - Test section 3; (1) heated- shield thermocouple, (2) multiple- shield thermo- 
couple, (3) Bureau of Standards gold-shield thermocouple, (4) pitot-static tube. 
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Figure 7. - Heated-shield component structure and assembly. (Drawn four times actual size.) 


Weight flow at station 2, W 2 , lb/sec 
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Figure 8. - Variation of weight flow with operating temperature. 


Mach number 
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Reference gas temperature, °R 

Figure 9. - Variation of Mach number with reference gas temperature under 

varying weight-flow conditions. 
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Indicated temperature 
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Reference gas temperature, °R 

Figure 10. - Variation of test-section-wall temperature from reference gas 

temperature for all conditions. 


1222 


1222 



p 

Radius, r , sq In. 

(a) Reference gas temperature, 550° R. 


Figure 11. - Variation of pressure profile at various test-section locations from pipe 
center line along diameter of pipe at constant reference gas temperature and constant 
weight flow of 6.3 pounds per second. 
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(b) Reference gas temperature, 1500° R. 


Figure 11. - Continued. Variation of pressure profile at various test-section 
locations from pipe center line along diameter of pipe at constant reference 
gas temperature and constant weight flow of 6.3 pounds per second. 
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Radius, r 2 , sq in. 

(c) Reference gas temperature, 2300° R. 

Figure 11. - Concluded. Variation of pressure profile at various test-section 
locations from pipe center line along diameter of pipe at constant reference 
gas temperature and constant weight flow of 6.3 pounds per second. 
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Figure 12. - Variation of C determined from pressure -profile surveys under varying gas-temperature operation. 
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Figure 13. - Evaluation of C by weight-flow considerations under varying gas-temperature operation. 
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Figure 14. - Comparison between weight-flow evaluation method and pre ssure -prof ile survey method for deter- 
mination of ratio f under varying gas -temperature operations. 
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Figure 15. - Variations of dimensionless coefficient N for all 
conditions of gas flow and reference gas temperature. 
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Figure 16, - Comparison of gas -temperature profile in 
test section 1 measured by thermocouple rake to gas 
temperature indicated by reference system. 
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(a) Weight flow, 4.6 pounds per second. 

Figure 17. - Variation of gas-temperature Indications from conven- 
tional thermocouples with gas -temperature indications from reference 
system. 
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Reference temperature, °R 

(b) Weight flow, 5,3 pounds per second. 

Figure 17. - Continued. Variation of gas-temperature indications from con- 
ventional thermocouples with gas -temperature Indications from reference 
system . 
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(c) Weight flow, 6.6 pounds per second. 

Figure 17. - Concluded. Variation of gas -temperature indications from con- 
ventional thermocouples with gas-temperature indications from reference 
system. 
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(b) Weight flow, 6.6 pounds per second. 


Figure 18. - Percentage deviations of conventional thermocouple 
temperature indications from pressure-sensitive system (reference 
temperature ) . 
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Figure 19. - Variation of gas -temperature indications for varying 
thermocouple -bead diameters at constant reference temperature. 


Indicated gas temperature 


50 


NACA TN 2043 






























































































Thermocouple 







/ 







O Multiple shiel 
p Average of the 
couple survey 
Q 3/8-in. diara 
Reference 

d 

irmo- 
r rake 






/ 















































/ 


^ 

n 

















J 

D U 

qQ 






1 ' 









o 

r~i 

[ 

£ 

] 















0 

-1 

O 









1 






C 

c 


















c 












> 




o 

a 
































) 

































l 1 


2800 


2600 


2200 


2000 


1800 


1600 


1400 


1600 1800 2000 2200 

Reference temperature, 


2400 


2600 


2800 


Figure 20. - Temperature deviation of conventional thermocouples from reference 
system with wall cooling. Weight flow, 6.6 pounds per second; wall temperature, 
710° R. 
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(a) Weight flow, 4,6 pounds per second. 

Figure 21. - Percentage deviations of conventional-thermocouple temperature 
indications from reference system with water cooling on walls. Wall tempera- 
ture, 710° R. 
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(b) Weight flow, 6.6 pounds per second. 

Figure 21. - Concluded. Percentage deviations of conventional-thermocouple 
temperature indications from reference system with water cooling on walls. 
Wall temperature. 710° R. 
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